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ABSTRACT 

Various models for electromagnetic emissions correlated with the gravitational wave signals expected to be 
detectable by the current and planned gravitational wave detectors are studied. The position error on the location 
of a gravitational wave source is estimated, and is used to show that it could be possible to observe the electro- 
magnetic counterparts to neutron star-neutron star or neutron star-black hole binary coalescences detected with 
the Advanced LIGO and the Virgo detectors. 

Subject headings: binaries: close — gravitational waves — techniques: miscellaneous 



1. INTRODUCTION 

A number of large laser interferometric detectors of gravi- 
tational waves (GW), developed by the LIGO project in North 
America, the GEO and the Virgo projects in Europe, and the 
TAMA300 project in Asia, are rapidly approaching their sensi- 
tivity goals. Months long data taking runs with three or more 
detectors are planned or have been completed in 2003. In ad- 
dition, a significant upgrade of the LIGO project interferome- 
ters which is planned for the end of the decade, the Advanced 
LIGO system, will lead to a ten-fold improvement in sensitivity 
in these detectors. Operating these GW detectors in a tightly 
coupled network has the advantage of reducing the likelihood 
of false detections, or equivalently of achieving a better distance 
reach for a given false alarm rate. The network can also be used 
to infer the position of the source. I argue in this paper that 
the precision of the position estimation and the electromagnetic 
(EM) fluxes expected from the most easily detectable sources 
of GW should be marginally sufficient to allow the observation 
of EM counterparts to GW events. 

Current estimates suggest that the most likely GW signal to 
be observed is the "chirp" from the in-spiral of two compact 
objects (neutron stars [NS], or black holes [BH] 1 ) in a close 
binary. A single (Advanced) LIGO detector should achieve 
a detection rate for NS-NS and NS-BH compact binaries of 
2 x 10" 3 -0.7 yr" 1 (10- 10 3 yr" 1 ) [Belczynski, Kalogera & 
Bulik (2002)]. These ranges correspond to a number of bi- 
nary formation models; the "Standard model" of (Belczynski, 
Kalogera & Bulik 2002) gives a NS-NS (NS-BH) coalescence 
rate of 50 Myr" 1 (10 Myr" 1 ) in the Milky Way. For NS-NS 
binaries, the signal-to-noise ratio (SNR) in a single LIGO de- 
tector (Advanced LIGO detector), assuming optimal signal pro- 
cessing, will be 10 for an optimally oriented source at 25 Mpc 
(425 Mpc), and will scale inversely with the source distance 
(Finn & Chernoff 1993). Averaging over source position and 
orientation reduces the distance for fixed SNR by a factor of 
~ 5/2. However, a network of three similar, independent de- 
tectors would achieve the same SNR at a distance larger by a 
factor of ~ v3. 

Only limited attention has been given to the important prob- 
lem of understanding the EM emissions of compact binary 
mergers. I review in section 2 a few mechanisms which may 
be important for such emissions. I then present in section 3 

1 Wherever it is relevant, the mass of NS is assumed to be 1 AMq , and the mass 
of BH, 10M Q . 



an overview of the techniques which will be used to analyze 
the GW data in order to discover and locate coalescing com- 
pact binaries. A comparison of the localization performances 
of these techniques and of the observational capabilities of EM 
detectors is presented in section 4, and is used to estimate the 
observation rate of counterparts to compact binary mergers, for 
various networks of interferometers. 

2. ELECTROMAGNETIC SIGNALS 

I review below three different mechanisms that may be im- 
portant in generating EM emission when a compact binary co- 
alesces: magnetospheric interactions, the radioactive decay of 
ejected material, and relativistic blast waves. These models are 
rather crude in their predictions, and it is not clear which, if any, 
may accurately describe actual EM counterparts. They should 
be used below as order-of-magnitude estimates of the EM sig- 
nal. Alternatively, the observational prospects quoted below 
can be interpreted as prospects for setting upper limits on these 
models from the EM observation of GW events. 

When a NS orbits a strong magnetic field companion, an 
electric field is induced in the orbiting star, leading to par- 
ticle acceleration in the form of a stellar wind, and coher- 
ent EM radiation as in normal pulsars. (Hansen & Lyu- 
tikov 2001) predict that this results in a radio burst pre- 
cursor occurring seconds before the GW burst from the 
merger and during the in-spiral, with flux at 400 MHz of 
F ~ 2.1 mJy(r/100 Mpc)" 2 (B/10 15 G) 2 / 3 , for r the distance 
to the source, and B its magnetic field intensity. In addition, 
a large amount of energy is extracted from the orbital mo- 
tion and released in the magnetosphere as Alfven waves and 
a pair plasma. This energy could drive a relativistic wind of 
pairs and photons, which would become optically thin after 
some expansion (still before the binary merger), and yield 
a X-ray thermal emission lasting a few seconds with flux 
F~3x 10" 9 erg cm" 2 s^r/lOO Mpc)" 2 (B/10 15 G) 2 , and tem- 
perature increasing from 10 keV to 100 keV during the burst, 
according to the same authors. The fact that the magneto- 
spheric interaction model leads to bursts that occur before the 
end of the GW signal is problematic for the X-ray burst, but, 
as pointed out by (Palmer 1993), interstellar dispersion could 
delay the radio burst enough to allow its observation. Alter- 
natively, if X-ray or radio instruments monitoring most of the 
sky are available, archival searches for coincidences with GW 
signals could be possible. 

In the second scenario, it is the radioactive decay of the neu- 
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tron rich nuclei of the decompressed NS matter ejected dur- 
ing the merger that produces the energy required to power an 
EM signal. Numerical simulations (Rasio 1999) suggest that 
a mass shedding instability ejects ~ 10% of the NS mass af- 
ter one or two orbits following first contact. (Li & Paczyriski 
1998) consider a simple model where some of this matter ex- 
pands in a spherical envelope of mass M, which is heated by 
the radioactive decay of the ejecta. If a fraction / of the en- 
velope mass decays and is converted to heat, they predict a 
peak luminosity L ~ 2 x 10 44 ergs s- 1 (//10" 3 )(M/0.01 Mg)" 1 / 8 
and an effective temperature at peak luminosity around 3 x 
10 4 K(M/0.01 MeO^f/lO" 3 ) 1 / 4 , so that most of the EM ra- 
diation will be emitted as soft UV. The optical remnant should 
decay on a timescale of <~ 1 day. 

Finally, a third possibility is that the merger of a compact bi- 
nary generates a relativistic blast wave, as it has been argued 
to explain long wavelength counterparts to gamma-ray bursts 
[e.g., van Paradijs, Kouveliotou & Wijers (2000)]. The black 
hole accretion disk model (Meszaros & Rees 1997) predicts 
that NS-NS or BH-NS coalescences result in an intermediate 
state where an accretion disk forms around a remnant BH, and 
that this accretion disk is responsible for the jet that powers the 
gamma-ray burst. The disk is accreted on viscous timescales, 
so that this model is likely to be useful only for short gamma- 
ray bursts [duration < 1 s, Narayan, Piran & Kumar (2001)]. 
Gamma-ray bursts, when seen on-axis if they are non-isotropic, 
are visible at cosmological distances, i.e. much farther out than 
where the GW signal from binary coalescences will be visi- 
ble; for such bright EM sources, simple time coincidences with 
GW detections should be sufficient to identify the EM counter- 
part (Kobayashi & Mezaros 2002). However, depending on 
the currently uncertain amount of beaming in short gamma- 
ray bursts, a possibly large fraction of GW events might only 
have "orphan" afterglows with no detectable gamma-ray emis- 
sion (Kobayashi & Mezaros 2002): the afterglow is produced 
by the decelerating outflow, so that the relativistic beaming of 
its radiation is decreasing with time, thus making the emission 
more and more isotropic, and thus more likely to be observed. It 
might also be that some compact binary mergers generate a fire- 
ball that powers an afterglow, but fails to generate gamma-ray 
bursts (Huang, Dai & Lu 2002), so the search for counterparts 
to GW bursts should not be limited to gamma-ray bursts. 

The detection of afterglows for short duration gamma-ray 
bursts appears to be much more difficult than for long duration 
gamma-ray burst. This could be due to larger error boxes for 
short bursts, or to intrinsically weaker afterglows. (Panaitescu, 
Kumar & Narayan 2001) argue for the latter hypothesis, scal- 
ing a model for long duration bursts to short duration ones by 
taking the energy of the short burst to be 5 x 10 51 /47r ergs sr" 1 , 
and the density of the surrounding medium to be low, 10~ 3 
cm" 3 , since binary mergers are likely to occur outside their par- 
ent galaxy. For a source at z = 1, the R band magnitude, 0.1 
day after the collapse, is then R ~ 24, and the 2-10 keV flux is 
Fi-wk ~3x 10" 13 erg cm" 2 s" 1 . The radio flux peaks <~ 3 days 
after the collapse, and is < 3 x 10~ 3 mJy 2 . 

Another source of GW with a well understood EM coun- 
terpart is the collapse of the core of massive stars in a super- 
nova explosion. According to recent estimates (Fryer, Holz & 

2 In extrapolating these numbers to smaller distances below, I assume a power 
law dependence of the optical and X-ray flux on the radiation frequency, F oc 
v~ p l 2 , for an electron energy index p ~ 2, and a synchrotron slope F oc u [ / 3 
for radio waves. 



Hughes 2001), a good channel for GW emission in stellar col- 
lapses is the excitation of a bar-mode instability in newborn, 
rapidly rotating neutron stars. From the data of (Fryer, Holz 
& Hughes (2001), Fig. 7), the collapse of a 15M Q star would 
produce a detection with SNR in a single Advanced LIGO de- 
tector of ~ 8(10Mpc/r)(/V c /100), where N c is the number of 
cycles over which the bar is assumed to remain coherent, and 
where optimal signal processing and orientation are assumed. 
The rate of core-collapse supernovae in the Galaxy is well- 
known, and lies between 0.007-0.02 per year. The fraction of 
all core-collapse supernovae that result in a neutron star with 
enough spin to develop the bar-mode instability is, however, 
largely unknown. Assuming this fraction to be one, averag- 
ing the response of the detector over the whole sky, and us- 
ing the extragalactic rate extrapolation method of (Kalogera et 
al. 2001), this gives an upper limit on the detection rate of 
0.3(N C / 100) 3 per year (SNR > 5). The mean maximum abso- 
lute B magnitude of Type II supernovae is -16.9, with standard 
deviation 1.4 (Miller & Branch 1990), so that the EM emission 
is detectable (B < 20, assuming no reddening) for any GW de- 
tection out to 10 Mpc. It should be noted that (Dahlen & Frans- 
son 1999) predict <~ 0.36 core collapse supernovae per square 
degree with R > 22, at any given time. For the farthest super- 
novae detectable using GW, supernovae unrelated with the GW 
signal might be observed optically within the source position 
error box. 

3. LOCALIZATION OF THE GRAVITATIONAL WAVE SOURCE 

The interest in the problem of efficiently detecting GW bursts 
in real data with a significant non-Gaussian noise component is 
currently driven by the availability of preliminary data from the 
LIGO and the GEO projects [e.g., Sylvestre (2002)]. Models 
of the data analysis algorithms that use Gaussian noise are very 
useful in getting good estimates of the performances of these 
algorithms, and should become more accurate descriptions of 
the real analysis as the quality of the data improves. With this 
caveat, it is currently well understood how the analysis of the 
data for GW signals of a precisely known form should be per- 
formed: given a set of parameters describing the signal, a tem- 
plate waveform is formed, and is correlated with the data. The 
full parameter space of the signal is explored with a finite num- 
ber of points chosen so that the mismatch between a true signal 
and its closest approximation leads to a reduction in SNR that is 
smaller than a few percents. This procedure, matched filtering, 
can be applied independently or coherently to geographically 
distributed detectors. In the former case, an estimate of the sig- 
nal parameters, including its arrival time, is obtained at every 
detector, and the source position is estimated by triangulation, 
using the relative phase of the signal in all detectors. In the 
latter case, the response to a template waveform is calculated 
coherently using the data from all detectors (i.e., by adding the 
log-likelihood from all interferometers for a single set of source 
parameters), and the position estimate is given by the point in 
parameter space with the largest response. The coherent ap- 
proach is more efficient in terms of detection efficacy and error 
boxe sizes (Pai, Dhurandhar & Bose 2001), but it might be 
computationally prohibitive to implement, at least in its most 
naive form (Pai, Bose & Dhurandhar 2002). 

In the results given below, I will consider the network of 
three interferometers consisting of the Virgo (V) instrument 
near Pisa (Italy), the LIGO instrument near Livingston (L), 
Louisiana (USA), and the 4 km LIGO instrument near Hanford 
(H), Washington (USA). In the case of triangulation, the form 
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of the position error, in terms of the solid angle Afl contain- 
ing the true source position 95% of the time, for experiments 
repeated on independent data but for the same true source posi- 
tion, is [Kip Thorne, cited in Gursel & Tinto (1989)]: 



Afltx 



c 2 At LV At 



LV^HV 



(1) 



A HLV \cos9\ ' 

where Af,j is the standard deviation on the time delay between 
interferometers i and j (the HL baseline can also be used if it re- 
duces the position error), Ahlv is the area of the triangle formed 
by the three interferometers, and 9 is the angle between the nor- 
mal to that triangle and the direction of propagation of the GW. 
Assuming Gaussian errors, (Afy) 2 = (Af,') 2 + (Af,) 2 , for Af; the 
standard error on the signal's arrival time estimation in detector 
i. The SNR and Af, are calculated for each interferometer for 
a "Newtonian chirp" for NS-NS and NS-BH coalescences as in 
(Finn & Chernoff 1993), but with up-to-date noise spectra 3 . I 
have calibrated Eq. (1) using numerical simulations, so that the 
proportionality constant multiplying the right-hand side of the 
equation is <~ 5.0. 

The position error from a detection of a NS-NS coalescence 
using the coherent approach is taken from (Pai, Dhurandhar & 
Bose 2001), and is given by 

A „ 2x3.7x lO^sr / 12 \ 2 

An = 1 m — > W 

|cos0| \Pn J 

where p^ is the network SNR, which is the sum in quadrature 
of the SNRs in all detectors. I include the extra factor of 2 to 
convert the quoted ler result to a 2a error box (95% probabil- 
ity coverage). With the initial network, the coherent method 
yields an error box the same size as triangulation for a source 
at 9 = 0, but a factor of ~ 3.6 farther (in good agreement with 
An oc SNR -2 ). Networks with detectors having characteristics 
similar to those of the Advanced LIGO detectors have not been 
studied by (Pai, Dhurandhar & Bose 2001), so I will interpolate 
their results by using Eq.2 with the network SNR computed for 
the Advanced LIGO noise characteristics, when detectors at all 
sites (including Virgo) are assumed to have noise performances 
close to the Advanced LIGO design. I will also consider the 
intermediate network consisting of two Advanced LIGO detec- 
tors and of the Initial Virgo detector. In that case, the error box 
will be an elongated ellipse, instead of being mostly circular. To 
account for this asymmetry, I use p 2 N = p v V ' p 2 H + p 2 L in Eq. (2), 
where the SNRs for the LIGO instruments are computed with 
Advanced LIGO noise power spectra, and the SNR for Virgo 
uses its Initial noise spectrum. The advantages of using the co- 
herent analysis are not as large for this network as for networks 
with similar instruments at all sites. 

4. OBSERVATION PROSPECTS FOR COUNTERPARTS 

Figure 1 presents the position error for a NS-NS binary 
merger as a function of its distance, the brightness of different 
R band models, and the rate of coalescence from the Standard 
model in a spherical volume extending out to that distance, us- 
ing the extragalactic rate extrapolation method of (Kalogera et 
al. 2001). The position error is a strong function of the posi- 
tion of the source on the sky, and the results of Fig. 1 are for a 
source located along the normal to the detectors' plane (9 = 0) 
and seen face-on. The right-hand side axis was scaled so that 

3 Initial LIGO: www.ligo.caltech.edu/~kent/ASIS_NM/noise_models.html; 
Virgo: www.virgo.infn.it/senscurve; Advanced LIGO: LIGO Document LIGO- 
M990288-A-M. 



the limiting magnitudes of two optical detectors correspond to 
their field-of- views. The first one is ROTSE-III [Smith et al. 
(2002), uppermost dotted line in fig. 1], a robotic telescope 
with limiting magnitude R <~ 18.5 for 1 minute of integration, a 
large field-of-view (3.4 deg 2 ), and which can respond in <~ 10 s 
to a trigger. The second one is the Suprime-Cam camera of the 
8.2 m Subaru Telescope (lowest dotted line), which has a large 
field-of-view (0.25 deg 2 ) for a limiting magnitude of R <~ 26 for 
a 10 minutes exposition 4 . 

In addition to this figure, Table 1 presents information for 
a full-sky coverage. The value quoted in the fourth column is 
the distance at which the position error for a given source and 
GW detection system is equal to the EM detector field-of-view 
(third column), assuming a binary seen face-on along the nor- 
mal to the detectors' plane. The fifth column presents the SNR 
in the EM detector for a source at that distance. For the ROTSE- 
III and Subaru telescopes, the radioactive decay model is used 
for the R band luminosity. For X-ray observations, I assume the 
use of the High Resolution Camera on-board the Chandra X-ray 
Observatory 5 , which has a 0.25 deg 2 field-of-view, with a lim- 
iting flux sensitivity of 6 x 10" 14 erg cm" 2 s" 1 for a 10 minutes 
integration time. The EM signal model is then the relativis- 
tic blast wave. Finally, I used the magnetospheric interaction 
model for radio observations, and assume the use of the Very 
Large Array telescope; in the P band, it has a 1 mJy RMS sen- 
sitivity for 10 minutes of integration, and a 5 deg 2 field-of-view 
(primary beam full width at half power of 2.5 deg) 6 . I assume 

4 www.naoj .org/Observing/Instruments/SCam 

5 hea-www.harvard.edu/HRC 

6 www.aoc.nrao.edu/vla/obstatus/vlas/vlas.html 
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FIG. 1. — The detectability of optical counterparts to the coalescence of a 
NS-NS binary, as a function of the luminosity distance of the source. The con- 
tinuous lines show the expected position error (left-hand side axis) using trian- 
gulation, and the dash-dotted lines show the error for a coherent analysis. The 
thin, medium, and thick lines are for the initial LIGO- Virgo network, the net- 
work with Advanced LIGO detectors, and the network with Advanced LIGO 
detectors and an "Advanced" Virgo detector with noise levels comparabale to 
those of Advanced LIGO detectors, respectively. The lines are terminated at 
the distance where the SNR from an optimally oriented source is below 5 in 
all detectors for triangulation, or the network SNR is below 5 for a coherent 
analysis. The lower dashed line shows the R apparent magnitude (right-hand 
side axis) expected at peak luminosity from the radioactive decay model. The 
upper dashed line shows the afterglow luminosity in the shock wave model, 0. 1 
day after the merger. The top axis shows the expected number of coalescences 
out to a certain distance. 



4 



Sylvestre, J. 



Table 1 

Observation prospects for counterparts 



source 


GW net a 


EM instrument 


distance 


EMSNR 


rate 








(Mpc) 




(yr" 1 ) 


NS-NS 


I 


VLA 


30 


30 


2 x 10" 3 


BH-NS 


I 


ROTSE-III 


40 


100 


3 x 10" 3 


NS-NS 


II 


ROTSE-III 


60 


3 x 10 4 


0.01 


NS-NS 


II 


VLA 


70 


4 


0.02 


BH-NS 


II 


ROTSE-III 


100 


900 


0.02 


NS-NS 


III 


Subaru 


100 


2x 10 4 


0.1 


NS-NS 


III 


Chandra 


100 


3 x 10 4 


0.1 


NS-NS 


III 


ROTSE-III 


400 


8 


7 


NS-NS 


III 


VLA 


500 


0.07 


10 


BH-NS 


III 


Subaru 


300 


4x 10 3 


0.2 


BH-NS 


III 


Chandra 


300 


6x 10 3 


0.2 


BH-NS 


III 


ROTSE-III 


900 


0.4 


10 



a The symbols I, II, and III refer to the initial network, the Advanced 
LIGO- Virgo network, and the network with Advanced LIGO detectors and 
an Advanced Virgo detector, respectively. In all cases, a coherent analysis 
is assumed. 



below in the rate estimates for radio observations that all NS bi- 
naries which coalesce contain at least one NS with a very strong 
magnetic field. Depending on the formation scenario of NS bi- 
naries, only a fraction of the binaries may actually contain a 
strong field NS, and therefore be candidates for radio emission 
by the magnetospheric interaction model. 

For every sky position, binary inclination angle, and GW po- 
larization angle, there is a maximum distance at which the posi- 
tion error box size equals the field-of-view of the EM detector. 
Integrating the binary coalescence rate over the whole sky and 
out to that maximal distance, and averaging over the inclination 
and polarization angles, gives the expected number of counter- 
parts that will be observable with that EM detector, per unit 
time. This rate is quoted in the last column of Table 1. It should 
be noted that the counterpart observation rates do not depend 
directly on the models of EM emission discussed previously. It 
is only assumed in computing these rates that the EM counter- 
part will be detectable out to the maximal distance where the 
error box size equals the field-of-view of the EM detector. For 
the models of EM emission under consideration, the fifth col- 
umn of Table 1 shows that this condition is nearly always met. 

If we happen to observe one NS-NS (NS-BH) merger with 
the initial network (by luck, or because the rate is larger by a 
factor of <~ 100 than the Standard model), it will have to be 
within 20 Mpc (40 Mpc) from us in order to allow its observa- 
tion in the optical wavelengths by a large field-of-view instru- 
ment like ROTSE-III, or within 30 Mpc for radio observations, 
and it will have to be near the normal of the detector's plane. 
With Advanced LIGO detectors, the predicted EM counterpart 
observation rate for NS-NS or BH-NS events is still only one 
every 25 years. The required EM detector's field-of-view to 
achieve a certain detection rate scales like the rate to the 2/3 
power 7 , so EM couterparts observation rates of one per year 
might be achieved with Advanced LIGO and Initial Virgo de- 

7 (detection rate) oc (distance reach) 3 , and (detector's field of view) ~ (error 
box size) oc (SNR)~ 2 oc (distance reach) 2 oc (rate) 2 / 3 



tectors if it were possible to effectively detect the EM source 
in a region of <~ 40 deg 2 . With a futuristic GW network where 
Virgo has a noise level comparable to Advanced LIGO instru- 
ments, tens of observations of counterparts could be made every 
year. Occasional deep observation of binary mergers might also 
be possible with 8-meter class optical telescopes and the Chan- 
dra X-ray observatory. Given the errors in the estimate of the 
compact binary merger rate, the predicted rates of counterpart 
observations could be an order of magnitude larger or smaller. 
Consequently, the observation of EM counterparts to compact 
binary mergers is improbable for the Initial network, possible 
with Advanced LIGO detectors, and likely with an Advanced 
detector in Europe. 

The optimal exploitation of the world-wide GW network will 
require the use of coherent analyses, and the development of 
advanced computational strategies might be necessary for this 
to be possible. Fully coherent analyses might be triggerred by 
less expensive incoherent ones, and therefore be ran only on a 
small subset of the data. In most cases, the EM counterparts 
might be rapidly dimming objects, so that it would be desirable 
to initiate the observation campaign within a few hours of the 
GW signal, at most. It is currently planned that GW data will 
be analyzed nearly in real-time, and the rapid transfer of data 
from multiple sites to a central processing facility has already 
been demonstrated (Marka, Mours & Williams 2002). The 
development of a world-wide warning system [like SNEWS, 
Scholberg (2000)] could also prove essential. In addition, some 
sources, like supernovae, might never have a GW waveform 
known with enough sophistication to perform matched filtering. 
It will then be necessary, in order to observe counterparts to 
these GW events, to develop robust source positioning methods 
for GW networks [Giirsel & Tinto (1989); Sylvestre (2003)], 
which do not require a signal model as precise as for matched 
filtering. 

The addition of other sensitive instruments in the GW net- 
work, especially in Asia or Oceania, would improve the like- 
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lihood of observing EM counterparts, not necessarily by dra- 
matically reducing the size of the best error boxes, but rather 
by improving the fraction of the sky where the error boxes are 
small enough to allow EM counterparts observations. Finally, 
it might be necessary to get beyond the rough EM models re- 
viewed in this communication in order to optimize the search 
strategies, and to correctly interpret future observations, or ab- 
sence thereof. Large regions of the sky will have to be searched 
in response to GW triggers, perhaps slightly beyond the capa- 
bilities of currently available observatories, especially with re- 
spect to the rejection of the astronomical background of vari- 
able objects. 

This work was supported by the National Science Founda- 
tion under cooperative agreement PHY-9210038 and the award 
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ratory document number LIGO-P020024-01-R. 
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